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We investigate the influence of the Coulomb interaction on the adiabatic pumping current through a quantum 
dot. Using nonequilibrium Green's functions techniques, we derive a general expression for the current based 
on the instantaneous Green's function of the dot. We apply this formula to study the dependence of the charge 
pumped per cycle on the time-dependent pumping potentials. Motivated by recent experiments, the possibility 
of charge quantization in the presence of a finite Coulomb repulsion energy is investigated. 
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I. INTRODUCTION 

The basic idea of electron pumping, put forward in the pi- 
oneer work of Thouless,— is to generate a DC current through 
a conductor in the absence of an applied bias voltage. This 
may be accomplished by applying time-dependent perturba- 
tions to the conductor. In electronic transport through meso- 
scopic conductors, the typical experimental time scale over 
which these external perturbations vary is large compared to 
the lifetime of the electron inside the conductor (dwell time). 
In that case, the pumping mechanism is called adiabatic. 

Adiabatic quantum pumping in mesoscopic noninteract- 
ing open quantum dots was investigated theoretically by 
Brouwer 2 by means of a scattering approach. Applying the 
emissivity theory introduced by Biittiker and co-workers, 3 
he demonstrated that the pumping current is proportional to 
the driving frequency and shows large mesoscopic fluctua- 
tions accounted by Random Matrix Theory. This scatter- 
ing approach has been employed to investigate several as- 
pects of adiabatic quantum pumping in noninteracting sys- 
tems, such as the role of discrete symmetries on the pumped 
charge, 4 the effects of inelastic scattering and decoherence, 5,6 
the role of noise and dissipation, 7 Andreev interference effects 
in the presence of superconducting leads, 8,9 as well as spin 
pumping . 10 ' 11 ' 12 - 1 3 Pumping phenomena in noninteracting sys- 
tems have also been investigated using alternative theoretical 
approaches, such as the formalism based on iterative solutions 
of time-dependent states 14 and the Keldysh formulation. 15 
Both approaches can be used beyond the adiabatic approxi- 
mation. 

Experimentally, the first implementation of an electron 
pump was due to Pothier et al. when charge was quantized 
due to Coulomb blockade (CB) effects. 16 Adiabatic phase- 



coherent charge pumping, though not quantized, was observed 
in open semiconductor quantum dots 17 and in carbon nan- 
otube quantum dots. 1819 Quantized charge pumping was re- 
cently observed in AlGaAs/GaAs nanowires using a single- 
parameter modulation, 20 a result with potential applications 
to metrology. An experimental realization of a quantum spin 
pump has also been implemented^ 

Pumping through interacting systems, where the scatter- 
ing approach does not apply, has been much less studied so 
far. Using the slave-boson mean-field approximation, Aono 
investigated the spin-charge separation of adiabatic currents 
in the Kondo regime. 22 The behavior of the pumping cur- 
rent through a quantum dot in the Kondo regime was stud- 
ied both for adiabatic 23 and nonadiabatic systems 24 using the 
Keldysh formalism. Quantum pumping was investigated both 
in the CB regim e) 25 ' 26 as well as for almost open quantum 
dots* 2 - 7 - The nonequilibrium Green's functions technique has 
been employed to investigate adiabatic pumping through in- 
teracting quantum dots in infinite U systemsi^^ 9 . The role 
of the Coulomb interaction in the adiabatic pumping current 
has also been investigated in the limit of weak tunneling and 
infinite-?/ using diagrammatic techniques. 30 The presence of 
electron-electron interactions was shown to improve charge 
quantization in one-dimensional disordered wires under cer- 
tain circumstances. 31 The effects of the coupling of the quan- 
tum dot to bosonic environments and its implications to charge 
quantization were analyzed in Ref . [32|. The interplay of nona- 
diabacity and interaction effects on the pumping current were 
also recently reported* 3 ^ 4 . 

In the present paper we investigate adiabatic charge pump- 
ing through interacting quantum dots in the CB regime for 
temperatures much higher than the Kondo temperature. We 
consider quantum dots with a single level subjected to a fi- 
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nite Coulomb repulsion U in the case of double occupancy. 
We investigate the time dependence of the pumping current 
by keeping U finite, a scenario out of the domain of validity 
of the theory developed in Refs. |28|and|29|. This allows us to 
identify the relevant time scales controlling the current ampli- 
tude in realistic situations. We develop a general formalism, 
based on non-equilibrium Green's functions, to investigate the 
influence of the Coulomb interaction on the adiabatic pump- 
ing current. We discuss some applications and consequences 
of this formulation and evaluate several quantities of interest 
numerically for a range of parameters. Finally, the possibil- 
ity of charge quantization in the presence of a finite Coulomb 
repulsion is investigated. The study of charge quantization in 
the adiabatic regime is interesting by its own, and is also a nec- 
essary step towards the understanding of recent experiments 20 
dealing with non-adiabatic pumping. 

This paper is organized as follows. In Section|II]we present 
the model used to calculate the time-dependent current flow- 
ing through the quantum dot. Section [Til] is devoted to the 
explicit calculation of the relevant Green's functions. In Sec- 
tion [IV] we apply this calculation to derive an expression for 
the pumping current in the adiabatic approximation for sys- 
tems with finite U. The numerical evaluation of the current 
as well as a discussion of its consequences and implications 
is presented in Section [V] Finally, Section [VI] is devoted to a 
brief summary of our findings and concluding remarks. 



H. MODEL FOR TRANSPORT IN QUANTUM DOTS 

We consider a quantum dot (QD) with a single, isolated 
resonance in the Coulomb blockade regime, as schematically 
depicted in Fig.[T] The potential in the dot is controlled by a 
time-dependent gate voltage V g (i) such that the QD Hamilto- 
nian reads 



Hdat = E £ si t ) d t d s + Un^ni, 



(1) 



8=U 



where n s = d\d s is the number operator and d\ (d s ) is the 
creation (annihilation) operator for an electron with energy 
e s (t) = eos — i] eV g (t) and spin s in the QD. Here, e denotes 
the electron charge and rj is a lever arm factor for the gate 
voltage. Two single-channel leads are attached to the QD. It 
is assumed that electrons in the leads are noninteracting and 
obey the Hamiltonian 



-Hlead = 



E 



E E 

:=L,Rs=T,| 



Skas c kas C ka 



(2) 



where ct and c kas are, respectively the creation and annihi- 
lation operators for electrons with momentum k and spin s in 
the lead a. The QD is separated from the leads by tunneling 
barriers controlled by the lateral gates V\ and Vi (see Fig . [TJ . 
The coupling Hamiltonian reads 



lead — dot 



E 



• f £ 
kas L 



H.c. 



(3) 
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FIG. 1: (a) Schematic view of a two-contact quantum dot coupled to 
a time-dependent gate, (b) Sketch of the energy levels of the model 
described in the text. 



The tunneling matrix elements Vk a connect states in the leads 
to the resonant state in the dot and are assumed to be spin 
independent. The total Hamiltonian of our model is the sum 
of these three contributions, 



H = -ffload + Hdot + Hl ca d- 



dot- 



(4) 



The coupling between the states in the leads and those in 
the dot, combined with the dot charging energy, turns the 
time evolution of the system into a nontrivial many-body 
problem. As a result, we cannot apply a single-particle 
formalism to describe the transport through the system and 
the usual scattering-matrix formulation for pumping currents 2 
is inappropriate. To circumvent these difficulties, we em- 
ploy the Schwinger-Keldysh formalism and the equation-of- 
motion method 35 to calculate the current through an interact- 
ing quantum dot in the CB regime. 

Our starting point is the general expression for the time- 
dependent current in terms of the quantum dot Green's func- 
tion G S)S (t, t')m- 

J a {t) = ~lm{ f dt'V^{t')e^^l h V ka {t) 

n k,s J -°° 

x [f a (e kas )G r s Jt,t') + G< s (t,t')] }, (5) 



is the Fermi function 



where f a (E) = [ e ( B ^ Q )/fc B T + 
for the lead a maintained at a chemical potential /i Q and tem- 
perature T and ks is the Boltzmann constant. Throughout the 
text we consider pumping in the absence of an external bias, 
that is, /i = hl = £p. For convenience, we set £p = 0. 
The lesser, retarded, and advanced dot Green's functions are 
defined as 35 



G< s (t,t<) = -<4(f'K(*)>, 

g; s (m') = -le(t-t')({dM4(m, 

G a s Jt,t>) = l -6{t' -t)({d s (t),dl(t')}). (6) 

Now it remains to compute the Green's function G StS (t,t') 
which involves the quantum dot states. This is where the 
many-body aspects of the problem make their way into the 
pumping current. Section[III]is devoted to this issue. 
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III. CALCULATION OF G a 



which obeys the equation-of-motion 



The current in Eq. © is given in terms of the quantum 
dot Green's functions G r s s (t, t') and Gf s (t, t 1 ). To write ex- 
pressions for them, we start by calculating the time-ordered 
Green's function G s . s (t, t'), defined as^ 



(7) 



where T is the time-ordering operator. The equation-of- 
motion for G s .s is 



G s , s (t,t') = 5(t-t') + UG%(t,t') 



k , a 



In Eq. (O we have introduced the "contact" time-ordered 
Green's function 



G St kas(t, t') 



-H (T 



d s (t)c{ as (t') 



(9) 



(-ihJL - e kas ^ G s , kas (t,t') = V k * a (t')G s>s {t,t'), (10) 



as well as the second-order correlation function 



G?l{t,t') ee -UT[d s (t)n s (t)dt(t')]), (11) 



that involves four fermionic operators and is generated by the 
interaction term f/n^rif. The same interaction term leads to 
the appearance of even higher order correlation functions in 
the equation-of-motion for G^ , namely, 



ih—-e s {t)-U 
at 



G%{t,t') = 5{t-t'){n- s {t))+Y. [t& r 8L('>*0] +E h* r 2;L(M') - V^ltL(M') 

ka. ko> 



, (12) 



where the occupation number is defined as 

(n s (t)) = (dl(t)d s (t))=iHG< s (t,t) (13) 
and we have introduced three lead-dot correlation functions, 36 



ri3« i (< ) t , ) = -i(r[ Cfa „(t)n ff (t)dt(f)]) 



and 



T 



T 



4 aS (tK(t)d s (t)dl(t') 



c kaS (t)dt(t)d s (t)dW) 



(14) 



(15) 



(16) 



At this level, one can verify that the equations-of-motion do 
not close. Going to the next level, one obtains new (higher or- 
der) correlation functions and even more complicated expres- 
sions. To solve this problem, we shall recur to an approximate 
scheme, namely the mean-field approximation. 



A. Formal solution of the equations-of-motion within the 
Hartree approximation 

We now focus on the Coulomb blockade regime and ne- 
glect spin correlations in the leads. That is, we assume that 



the Kondo temperature,- T K ~ Uy/T/2U exp (-7r|e s |(£ s + 
U) /2UT) is very low, T K < T. As usual, T stands for the 
quantum dots resonance linewidth which will be precisely de- 
fined in Sec. [IV] Hence, with respect to Kondo correlations, 
we are in the high-temperature regime and the mean-field ap- 
proximation is expected to be valid. Within this approxima- 
tion, one can write the r( 2 )'s as 



rffifo*') = <»*(*)) G kas>s (t,t') 



and 



^(2)mf 

2;fco 



(*,0 = rgE(*.0 = o. 



(17) 



(18) 



It has been shown that Kondo correlations are still absent 
in the next order of the equations-of-motion hierarchical 
The latter dresses the Green's functions self- 



truncation 



36.3y 



energies with higher order terms in V that include, for in- 
stance, cotunneling processes. As long as e s is of the order 
of ksT, we have verified that these contributions give only 
small corrections to the Hartree mean-field approximation. 36 
Thus, we write 



ih^--e s (t)-U 
at 



p(2)mf 
^ ss,s 



(t,0 



K(t)) 



, (19) 
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where the occupation number (n s (t)) has to be determined 
self-consistently for all times. Equations ([8]), ( [Tol l, and ( fT9] l 
form a closed set of equations-of-motion that determines the 
time-ordered Green's function G s s . Using analytical contin- 
uation and the Langreth rule s 36 ' 40 we can then find the Green's 
functions G r s s and Gf s that appear in the expressions for the 
current, Eq. (|5). For convenience, let us define two auxil- 
iary time-ordered Green's functions g s and g^ 7 that obey the 
equations-of-motions 



-0 



g s (t,t') = 6(t-t') (20) 



and 



ih—-e.(t)-u 

at 



g u s {t,t')=S{t-t'), (21) 



respectively. By analytical continuation into the complex 
plane, we can rewrite Eq. ( fT9] > as 



G ( £f(T,T') = gV(T,T')(n,(T')) 
+ Y,J dr lffs c/ (T,T 1 )(n s -(T 1 ))V7 Q (T 1 )G feQS , s (r 1 ,r') 



(22) 



k, a 



The equation for Gkas,s (t~i , t') can also be obtained in a sim- 
ilar manner. Using Eq. ( fTUl ), the equation-of-motion for the 
time-ordered Green's function for free electrons in the leads, 
namely, 



or 



g kas (t,t') = 5(t-t'), (23) 



and the rules of analytical continuation, we conclude that 
the contour-ordered Green's function G s .k as (T, r') obeys the 
equation 

G Sy kas (r, t') = J dn G ss (r, n ) U fc * Q (n ) 5 fcct;5 (n , r'), 

(24) 

while its counterpart is given by 

Gfc as , s (r,T') = J dn g kas {T,Ti)Vk a (ri) G s>s (ti,t'). 

(25) 

In all these cases the integration paths run over the Keldysh 
contour discussed in Refs.l36landl4ll 

Now the equations-of-motions close since both G ( - 2 - )mf and 
Gkas,s are expressed in terms of G s . s and free Green's func- 
tions. By introducing the renormalized single-electron resol- 
vent 

g s (T,r')=g s (T,T') + U(ri s (T')) J dng s (T,n) g v s '(n,r'), 

(26) 

we write, after a little algebra, a Dyson-like equation for G StS , 



G S)S (t,t') = £f a (r,r / ) + J dri J dT 2 g s (T,T 1 ) 

xZ ss ( Ti ,t 2 )G s , s (t 2 ,t'), (27) 



with the self-energy defined as 

E ss (r,r') =J2 V ka^)9kas(T - T')V ka (r')- (28) 



The rather peculiar structure of our solution is noteworthy. 
The auxiliary Green's function g s , Eq. 1261 . is not a free prop- 
agator since it contains a term involving (n s ) that arises from 
the mean-field approximation and has to be calculated self 
consistently. The self energy carries information about the 
coupling to the leads and can be calculated independently of 
the state of the dot. Hence it does not contain information 
about the many-body character of the problem. 

In Section llVl we shall specialize the calculation to the adi- 
abatic regime, first by explicitly obtaining an expression for 
the Green's functions involved in Eqs. d26l l and (ffTJl and then 
by evaluating the current, Eq. (|5). 



IV. ELECTRONIC TRANSPORT IN THE ADIABATIC 
APPROXIMATION 

The two important time scales in the problem of charge 
pumping through non-interacting quantum dots are the mean 
dwell time of an electron inside the dot (lifetime of the reso- 
nant state), to, and the inverse of the characteristic pumping 
frequency, r pump = 27r/w pump . In typical experimental se- 
tups, the pumping frequency u> pump lies in the range between 
10 MHz to 1 GHz^i For uj puiap /2n = 100 MHz, one has 
Tpump = 10 ns. The mean dwell time is given by the inverse of 
the resonance width T. To estimate it, let us first recall that the 
dot single-particle mean level spacing is A = 2irh 2 / (Am*), 
where A is the dot effective area and m* = 0.067m e for 
GaAs. We obtain A w 7.6 /ieV(/im) 2 /.A, where A is given 
in square microns. For the Coulomb blockade regime, typ- 
ical resonance widths are T — 0.01 — 0.1 A. As a result, 
Td = fi/r ~ 0.8 — 8 ns(fim) 2 j A for most devices. For A 
much smaller than 1 (/im) 2 , we find that r pump S> td. In this 
case we can safely employ the so-called adiabatic approxi- 
mation, which precisely relies on the fact that the time scale 
over which the system parameters vary is large compared to 
the lifetime of the electron in the dot. 



A. Adiabatic approximation for the Green's functions 

A convenient way to separate slow and fast times scales is 
to reparametrize the Green's functions as 



G(t,f)-G(t-f,£±^ 



(29) 



that is, the time variables are replaced by a (fast) time differ- 
ence St = t — t' and a slow mean time t = (t + t')/2. We 
implement the adiabatic approximation to lowest order by ex- 
panding the Green's functions up to linear order in the slow 
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variables, namely, 



(t-f, t±^j « G(t-t',t) 



t'-t\ dG . 



(30) 



t=t 



In what follows we formally write 

G(t-t',t) = G (0) (t-i',i") + G (1) (t-t',i), (31) 



where the zeroth order refers to equilibrium quantities, while 
the adiabatic contributions, linear in the slow time variable 
(and in our case proportional to the pumping frequency), are 
collected in the first-order correction. The accuracy of our 
approximation can be tested by inspecting higher-order terms. 
We will return to this issue in Sec. [V] when we present our 
results. 

Let us now describe how the approximate scheme works. 
Using the mean-time parametrization, we write Eq. d26i i as 



g s (t-t',i)=g s (t-t',i) + U(ris(i}) J dt l9a ft - £ U - ? , 



(32) 



Expanding g s in the slow variables as in Eq. d30b and taking the Fourier transform with respect to the fast variable, namely, 

g(uj, t) = d(t - t') g(t - t', t) exp[«w(t - <')], we obtain 



with 



and 



g.(w,*)=ffi 0) (a;,t)+g«(w,t) ) 
;(°) ( W , t ) = ff <°> { U ,i) + U( 4 0) (t)\ 9 W (W) f ) 9 f (0) ( w , * ) 



9t 



<9w 



<9t duj 



+ (4 0) (t)) 



<9w <9t 



(33) 



(34) 



(35) 



where 



<n s -(t)) = (4 0) (t)) + (4 1) (t) 



(36) 



is introduced following the same principle as the one de- 
scribed after Eq. ( f30b . 

Equation d35l > is further simplified by the fact that the lowest 
order corrections to terms involving g a and gY^ vanish for 
the retarded component. To demonstrate this, let us consider 
the retarded component 



gl s {t-t') = -- e(t-t')exp 



dtie s (ti) 



(37) 



Expanding e s (ti) around the mean-time t = (t + t')/2, 
namely, e s (<i) = e s (t) + e s (t)(ti — t) we obtain 



dt 1 e a (t 1 ) = e s (i)t + 0(e), 



(38) 



S othat^ 1)r(o) =^ (1)r(a) 



0.— This simplification shows 



the advantage of the mean-time parametrization, Eq. d29l ), 
with respect to other parameterizations, such as the one cho- 
sen in Ref . G8L 



After these simplifications, we obtain for the advanced and 
retarded components 

#^(u;,t) = #^)(o; ) t) + (4 0) (t)) 

X gl°^ a Xu,t)U& i0Ma Hu,t) (39) 



and 



3?™M = („W(t)) ai°)-W( Wf i)^(W)( Wl i) 



is 



x 3. 



a(4 0) (*)) a 



U 



(40) 



For the lesser components, we employ the fluctuation- 
dissipation theorem to write 



#<(w,t) = /(w)[#>,tl-#>,t) 

and apply the Langreth rules to Eq. (l35l l to obtain 



(41) 
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jU<(u,,t) = C/(n«(f))/H [g^(co,t)g^(u;,t) -0<°>> S ^ (O) >,<) 



dt duo 



{/(a;) [# <>, t)g" W i) - W >, i)] } 



C/(0)a 



(0)r 



MO) 



>,«)]} 



Here /(w) = [exp(fiw/fc s T) + l]" 1 . 

We proceed in the same way to obtain an expression for G s . s . The result is 



with 



and 



G., a (w,t) = G$(w,i) + GW(w,t), 
G« ( Wj t) = (w, t) + s« (a;, i)E s ( W , t)G<°) ( w , t) + < 0) («, t)E. (w, t)G&) ( Wj f) 



,(o) 



ffW(w,i)S fl (a;,i) 



<9G. 



Co) 



0i 



-(w,t) 



,(0) 



(42) 



(43) 



(44) 



(45) 



In Eq. ( 1451 ) we have introduced 

5W(«,t) = X;[v r te(*)^a(t)-H., 



9w 



(46) 



In what follows we use the wide-band approximation, where 
E(w, t) — > E(i), in which case the above equations are sim- 
plified further. 

From Eqs. ( l44b and ((45), we obtain G r and G < , which 
are needed to calculate J a , Eq. (O, in the adiabatic ap- 
proximation for the Coulomb blockade regime. Since the 
zeroth order terms are essentially equilibrium quantities, 
we are allowed to use the fluctuation-dissipation theorem 
to compute G (0)< without much effort: G (0)< (w,t) = 
-2i/(w)Im [G(°>(u,t)]. For G^< this is no longer pos- 
sible and we have to use the Langreth rules. The resulting 
expressions are rather long and will be omitted here. 

The occupation numbers (n^) and (n^) that appear in 
Eqs. (l44l and (|43T > are calculated self consistently using 



(*)) = 



where i = or 1 . In the absence of an external magnetic field, 



which is the case considered here, ( n$ 



Ai) 



For later convenience, we assume the couplings Vka to be 
energy independent and use the flat and wide band approxi- 
mation to define 

r a (e,i) = 2ir\V a {e,t)\ 2 p a = 2ir \V a (t)\ 2 p a =T a {t), 

(48) 

with p a denoting the density of states in the lead a. We also 
introduce 



(49) 



as the total decay width. As we discuss next, the current in 
Eq. (0 is easily cast in terms of these quantities. 



B. Current in the adiabatic approximation 

To evaluate the time integral in the general expression for 
the current, we proceed as in Eq. ( f30b and expand all terms 
in the integrand to linear order in the slow variables. The re- 
sulting expression for the pumped current depends explicitly 
on G < (lu, t) and G r (ui, t). Since G < is related to occupations 
(and hence to fluctuations) and G r to dissipation, as shown 
by standard linear response theory, it is natural to break the 
current into two parts, 

J a (t) = J«(t) + J* is (t), (50) 
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where the fluctuation term is 



V. RESULTS AND DISCUSSIONS 



s 

= -| r *(*)E <"•(*)> 

s 

while the dissipation term is given by 



r Q (t) 



2^ 



G< (w,t) 



(51) 



2e 



r Q (t)G;>,*) 

fO(M). (52) 



Now we are ready to use the adiabatic expansion for the 
Green's function, G S)S = G^l + Gi]l, and to identify the 
zeroth and the first-order contributions to the pumped cur- 
rent, j(°) and JW, respectively. It can be shown that the 
zeroth order current vanishes, as expected by the fluctuation- 
dissipation theorem. 

The first-order contribution to the current due to fluctuation 
is given by 



ra(t)£(»s (*))> 



(53) 



while the first-order dissipation term is given by 

J« dis (t) = J<>) dis (t) + 4 lb)dis (i), (54) 

where 

(55) 



and 



4 lb)dis (t) 



_d r 

di 



(56) 



The reason for breaking the dissipation term into two con- 
tributions is that Ja lb ' dls (t) is a total derivative in time. In- 
tegrated over a pumping period, this current term does not 
contribute to the pumped charge. This provides a good check 
for the numerical calculations presented in Sec. |V] We also 
successfully verified that our analytical expressions yield the 
same results as other pumping formulations 2 ' 28 in the U — > 
limit. 

Equations (1501 . (1531 . and d55l l constitute the principal re- 
sults of this paper. In the following, we will use these ex- 
pressions to investigate the role of interactions on the pumped 
current. Specifically, we will study how interactions affect the 
dependence of the pumped current on U, temperature, and the 
phase difference between the pumping perturbations. 



In this Section we compute numerically the pumping cur- 
rent, Eq. (1501 . and investigate the dependence of the magni- 
tude of the leading contribution to the total charge pumped 
per cycle, 



Q 



(57) 



on several model parameters. In particular, we discuss in 
which conditions the pumped charge can be quantized to its 
maximum value, |e|. To accomplish this goal, we consider the 
following parametrization for tunnel couplings: 



T a (t) = T 0iCC + AT a cos (fit + <f> 



(58) 



where a — R,L and Tq a and Ar a are real constants. We 
also assume that the quantum dot resonance energy varies in 
time as 



r(t) = e + £ i cos(Qi). 



(59) 



Notice that since e s = e s , we have dropped the spin index. 
In the following, all parameters are chosen to ensure that the 
system is clearly in Coulomb blockade regime, T <C U. Typ- 
ically, we take r , a /U = T /U = 0.1 and AT a = AT = 
0.05 U in our numerical calculations. 

As already stressed, the analysis is restricted to the first- 
order adiabatic correction. Hence, since the current is linear 
in ft, the charge pumped per cycle does not depend on the 
pumping rate. The accuracy of this approximation depends on 
the magnitude of the second-order corrections. Intuitively, the 
adiabatic approximation becomes more accurate as the ratio 
Ml /Tq becomes smaller. A closer analysis of the time deriva- 
tives of the Green's functions induced by the adiabatic expan- 
sion reveals that the dimensionless parameter controlling the 
adiabaticity is rather £ = max{/ifi/To, MIsi/Tq}, Albeit the 
fact that the results presented here are always valid for a suf- 
ficiently slow pumping, such that ( <C 1, there is no simple 
way to estimate the accuracy of the approximation for a given 
pumping rate ft. To be quantitative, one has to evaluate the 
second-order correction within the adiabatic approximation, 
which is a quite daunting task. Instead, we did a rough esti- 
mate of these higher-order contributions by studying a single 
representative term that appears in the second-order Green's 
function. We found that it scaled with £ as predicted, up to a 
numerical factor of order one. 

Figure|2]displays the result of the self-consistent calculation 
of the zeroth order occupation (n.^ ), Eq. d47| >, as function 
of the position of the resonance e for three temperature val- 
ues. Knowledge of (n^ ) is crucial for computing the various 
terms that enter in the calculation of the pumping current. As 
expected, the occupation of the quantum dot increases when- 
ever the position of any of its two levels, e and e + U, co- 
incides with Fermi level ep = 0, facilitating charge trans- 
port. For low temperatures, this is the dominant mechanism 
of transport, whereas for higher temperatures thermal fluctua- 
tions can also induce charge transfer through the quantum dot. 
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FIG. 2: Equilibrium quantum dot occupation number (rig ) as a 
function of the level position e for three values of the temperature: 
k B T/U = 0.01 (full black curve), k B T/U = 0.05 (dotted blue 
curve), and ksT/U — 0.1 (dashed red curve). Here To = 0.1 U 
and the Fermi energy is set to zero, ef = 0. 

This explains why the features in the curve become sharper as 
temperature decreases. 

The first-order correction to the quantum dot occupa- 
tion number {rig 1 ^), also calculated self consistently using 
Eq. d47l >. is shown in Fig. [3] as a function of time for several 
values of eq. It is important to emphasize that (n^ 1 ') is intrin- 
sically a time-dependent quantity and depends on the pump- 
ing parameters dynamics, in contrast to (rig). Notice that the 
magnitude of (n'P) is typically much smaller than (nf^). We 
observe that the maximum values of (rig 1 " 1 } occur for £o = ep. 
When the position of the level £o deviates significantly from 
Ef, charge pumping is attenuated and the magnitude of the 
current is smaller. 




0.0 0.5 1.0 1.5 2.0 



t x n [>r] 



FIG. 3: First-order correction to the quantum dot occupation number, 
(ri'P ) , as a function of time over a complete pumping cycle for three 
values of £o: £o/U — —0.075 (dotted blue curve), ea/U = (full 
black curve), and Eo/U = 0.075 (dashed red curve). Temperature is 
k B T/U = 0.01, = -4>r = vr/2, e 1 /U = 0.05, T /U = 0.1, 
and AF/U = 0.05. 

After computing (n^), the next step is to calculate the 



first-order correction to the time-dependent current Ja\t) 
given by the sum of the fluctuation term Ja fl (i), Eq. d53l >, 

and the dissipation terms ji la ^ dls (f ) and Ja b ^ dm (t), Eqs. ( l55l > 
and d56| i, respectively. A typical result is shown in Fig. [4] 
where we plot the frequency independent quantity /f2 as 
a function of time over a full pumping cycle. It is important 

to point out that the second dissipation term, Jas ls (t), does 
not contribute to the total charge pumped per cycle since it is 
proportional to a total time derivative. Consequently, its time 
integral over a complete pumping cycle must vanish, a result 
that has been confirmed numerically. The analysis of Fig. [4] 
reveals that these three current terms, as (n)p), exhibit max- 
ima precisely at the instants when the resonance energy level 
sit) crosses the Fermi energy. In the case of Fig. [4] where 
£o = 0, these maxima occur at t — 7r/2fi and t — 37r/2f2. 



0.4 - 




0.0 0.5 1.0 1.5 2.0 

FIG. 4: The three terms that contribute to the first-order correction 
to the pumping current as a function of time: (t) (dotted blue 
curve), Ja 13 (t) (dashed red curve), and Ja b ^ dls (t) (full black 
curve). Here we set so = and take the other model parameters 
as in Fig. [3] 

There is an intuitive interpretation for the role of the pump- 
ing parameters of our model, Tn i,(t) and e(t), that helps us to 
understand the time dependence observed above: In Eq. (T59b 
we fixed the phase offset of e(f) to zero. In this situation, 
for < t < Tpump/2 the resonance energy e decreases 
with time. As a consequence, during this half pumping pe- 
riod (n s ) increases with time, which corresponds to loading 
negative charge into the quantum dot. In this time interval, 
the sign of the pumping current depends on the phase dif- 
ference between (f>R and (p^. The situation is reversed for 
Tpump/2 < t < T pump . Figure|5]shows the three-dimensional 
plot of the charge pumped per cycle Q as a function of both 
(f>u and (j>L. Consistent with the reasoning presented above, 
having <f>L and <f>R in anti -phase favors larger values of \Q\. 
In particular, we find two maximum values of \Q\, one at 
4> L = 7r/2 and <j)R = 3tt/2, and the other at <pi = —ir/2 
and <f>n = 7r/2. The location of these maxima shows no de- 
pendence on any of the model parameters, provided e\ ^ 0. 
In this limit case, there are only two active pumping parame- 
ters, Tr and Tl, and the dependence of Q on the (f>R and 4>l 



9 



is the same as in the non-interacting case. 2 Since we are in- 
terested in maximizing \Q\, in the remaining of this paper we 
take E\ ^ and </>l = — (pR = t/2. 




FIG. 5: Three-dimensional graph of Q as a function of <j>L and <j!>_r. 
Temperature is k B T/U = 0.01 while e = 0, Ei/U = 0.05, 
To/U = 0.1, and AT/17 = 0.05. 

We are now ready to study the dependence of Q on V g (t), 
related to £q and e\, as well as on the dot-lead couplings, rep- 
resented in our model by To and AT. 

In Fig.[6]we show the charge pumped per cycle Q calculated 
as a function of eq. Charge pumping is enhanced whenever a 
quantum dot resonance, £o or £o + U, crosses the Fermi level, 
resulting in the two peaks of Fig. [6] 
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FIG. 6: Charge pumped per cycle as a function of the level position 

£ forei/(7 = 0.05, k B T/U = 0.01, (j> L = ~<t>R = tt/2, F /U = 
0.1, and AF/U = 0.05. The charge is measured in units of the 
electron charge e. 

Figure Q shows the dependence of Q on E\. We consider 
one of the situations of maximum pumping, namely, <f>L — 
—cf>R = ir/2 and £q = £f = 0. In this case, \Q\ increases 
monotonically with E\. We caution that once E\ exceeds To, 
it is necessary to check whether £ <C 1, so that the adiabatic 
approximation still holds. Hence, increasing e% might not be 
advantageous whenever it is necessary to reduce fl Figure [7] 
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FIG. 7: Charge pumped per cycle as a function of the resonance 
oscillation amplitude e x for e = 0, F /U = 0.1, AF/U = 0.05, 
kBT/U = 0.01, and 4>l = — <f>R = n/2. The charge is measured in 
units of the electron charge e. 

also shows that Q vanishes when E\= 0, as expected for a two- 
parameter adiabatic pump that occur for 4>l = — <f>R = tt/2. 1 ' 2 
We now address the dependence of Q on Ar and r . To be 
quantitative, we now also keep T 3> Tk for the sake of the 
validity of our approximation. To maximize pumping, we find 
that it is advantageous to decrease Tk by taking £o / rather 
than increasing T. As before, we consider 4>l = —<f>R — 
tt/2. Due to the time derivatives appearing in the Green's 
function expressions, several terms in Eqs. (|53l l and d55l l are 
proportional to AT. Indeed, we find that Q is roughly linear 
in AF for several values of E\ < Tq. Figure[8]shows Q versus 
To for three temperature values. Due to the fact that IcbTk < 
\JTU /2e~ 7 '/ 2 for e w To, our approximation scheme breaks 
down as T is increased and T/< reaches T. 



-0.10 




0.0 0.1 0.2 0.3 0.4 0.5 



r /u 



FIG. 8: Charge pumped per cycle as a function of To for differ- 
ent values of temperature: kBT/U = 0.05 (full black curve), 
k B T/U = 0.1 (dotted blue curve), and k B T/U = 0.2 (dashed 
red curve). Here e = T, ei/U = 0.05, AF/F = 1, and 

4>L — ~4>R — 7t/2. 

Figure|9]shows Q as a function of temperature for three val- 
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ues of the resonance energy Eq with ej kept fixed. The tem- 
peratures for which we observe the largest values of \Q\ scale 
with £o- We also find that by decreasing |eo| the maximum of 
\Q\ increases. Unfortunately, since our results are only valid 
for T ;» Tk, we cannot freely vary e . 




0.00 0.05 0.10 0.15 0.20 0.25 0.30 



FIG. 9: Charge pumped per cycle as a function of temperature for 
E /U = 0.075 (full black curve), s /U = 0.1 (dotted blue curve), 
and Eo/U — 0.15 (dashed red curve). For all curves Ei/U = 0.05, 
F /U = 0.1, AF/U = 0.1, and cj> L = -<f> R = tt/2. 

Finally, let us address the dependence of Q on the charging 
energy U. Our results are summarized in Fig.QI)] A large in- 
terval range for U is displayed to best illustrate the pumped 
charge dependence on this parameter. We observe that pump- 
ing is largely enhanced for small values of U . When U be- 
comes comparable to T the system departs from the Coulomb 
blockade regime. 
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FIG. 10: Charge pumped per cycle as a function of U/Fo for eo = 
To, ei = To/2, and k B T = Ar = r . 



VI. CONCLUSIONS 

In conclusion, we have investigated adiabatic charge pump- 
ing through quantum dots in the Coulomb blockade regime. 
We specifically studied the impact of Coulomb interaction 
on the pumping current amplitude for the finite-U Anderson 
model, in contrast to previous works that treated the infinite- 
U case M 

We have derived a general expression for the adiabatic 
pumping current that is proportional to the instantaneous 
Green's function of the dot. This formula was then applied 
to compute the time dependence of the total charge pumped 
per cycle through the dot. This allowed us to analyze several 
aspects of experimental relevance, such as the dependence of 
the pumped charge on temperature and on the phase difference 
between time-dependent perturbations. 

We find that, within the adiabatic regime, there is a large 
range of parameters that can be used to maximize the charge 
pumped per cycle. For this purpose, we find that it is advan- 
tageous to: (i) Tune the back gate voltage to pump with the 
QD in resonance with the Fermi energy in the leads; (ii) max- 
imize the pumping amplitude AT and, possibly, E\ as well; 
(iii) minimize temperature. 

We were not able to find a set of parameter values that gives 
one unit of charge e per pumping cycle within the parame- 
ter ranges allowed by our approximations. We do not discard 
such interesting possibility, but our investigations hint that it 
may only be possible for very particular pulse formats, not 
necessarily sinusoidal, and within a narrow parameter inter- 
val. The possibility of spin pumping and the consideration 
of the double-dot case are under investigation and will be re- 
ported soon. 
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